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ABSTRACT 

If a microlensing event is caused by a star, the event can exhibit change in color due to the light from the lens. 
In the previous and current lensing surveys, the color shift could not be used to constrain the lens population 
because the blended light responsible for the color shift is mostly attributed to nearby background stars rather 
than the lens. However, events to be observed in future space-based surveys do not suffer from blending and 
thus the color information can be used to constrain lenses. In this paper, we demonstrate the usefulness of 
future surveys in measuring color shifts. By conducting simulation of galactic lensing events based on the 
specification of a proposed space-based lensing survey, we estimate that the shift in the color of R-H will be 
measured at Scr level for ^ 12% of events that occur on source stars with apparent magnitudes brighter than 
J = 22.5. Color-shifted events tend to have high magnifications and the lenses will have brightnesses equivalent 
to those of source stars. The time scales of the color-shifted events tend to be longer than those without color 
shifts. From the mass distribution of lenses, we find that most of the color-shifted events will be produced by 
stellar lenses with spectral types down to mid M-type main sequence stars. 
Subject headings: gravitational lensing 
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1. INTRODUCTION 

One important characteristic of a microlensing event is that 
the color of the lensed star does not change during the lens- 
ing magnification. However, events observed in actual lens- 
ing surveys can exhibit color shift. This color shift occurs 
when the flux from the lensed source star is blended with the 
flux from other stars located in the seeing disk of the star. 
The lensed and blended stars generally have different colors. 
Then, color shift occurs as the source/blend relative bright- 
ness changes during lensing magnificati on. 

Kamionkowski (1995) and Buchalter et al.l (Il996h pointed 
out that if the object responsible for a lensing event is a star, 
it is possible to obtain information about the lens from the 
measured color shift. However, this idea could not be imple- 
mented for the actual analysis of lens population. The most 
important reason for this is that under the previous and current 
lensing surveys conducted by using ground-based telescopes, 
the blended flux is mostly attributed to nearby background 
stars rather than the lens. 

However, the situation will be different in future space- 
based lensing surveys, such as the Microlensing Planet Finder 
(MPF). The MPF is a proposed space mission to NASA's 
Discovery Program exclusive for microlensi ng and it will be 
equipped with a 1.1 m a perture telescope (iBennett & Rhid 
l2002l;lBennett et alj2004 . Such a survey is ideal for the mea- 
surement of the color shift in lensing events due to the follow- 
ing reasons. First, thanks to the high resolution provided from 
space observation, the survey does not suffer from blending 
caused by background stars. Then, one can safely attribute 
the observed color shift to the lens. Second, the fraction of 

'corresponding author 



color-shifted events would be higher in the space-based sur- 
vey because it will monitor much fainter stars than current sur- 
veys. As the monitored star becomes fainter, the relative flux 
of the lens to the source star increases and thus the event is 
more likely to produce large color shift. Third, with high pho- 
tometric precision, the survey can detect much smaller color 
shift than ground-based surveys. With the increased fraction 
of large color-shifted events and higher sensitivity to small 
color shift combined with nearly an order higher event rate, 
the future space-based lensing survey would be able to detect 
color shifts for a large number of events. 

In this paper, we demonstrate the usefulness of future 
space-based lensing surveys in measuring color shifts of lens- 
ing events. For this, we estimate the fraction of events with 
measurable color shifts among the events to be detected from 
the proposed MPF survey. We also investigate the properties 
of the color-shifted events by conducting detailed simulation 
of galactic lensing events under realistic observational condi- 
tions. The fraction of color- shif ted events was previously es- 
timated bv Buchalter et al.l ( 11996) . However, their work was 
intended as a template for general comparison with various 
observational assumptions and thus is not focused on a spe- 
cific survey. In addition, they did not provide information 
about the properties of color-shifted events. 

The paper is organized as follows. In § 2, we introduce for- 
malism used for the description of the color shift in general 
cases and under various limiting cases. In § 3, we describe 
the simulation conducted to estimate the fraction of color- 
shift events and their properties. In § 4, we present the result- 
ing distribution of color shift. We probe the properties of the 
color-shifted events by investigating the distributions of var- 
ious quantities related to the events including the source and 
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lens brightness, lens light fraction, magnification, lens mass, 
and event time scale. In § 5, we summarize and conclude. 

2. FORMALISM 

The difference between the colors of a lensing event mea- 
sured during and before lensing magnification is expressed as 

Fr/Fh 



A(/?-//) = -2.51og 



(1) 



^FRfi/Fufi 

where F represents the observed flux, the subscript '0' de- 
notes the measured quantity without lensing magnification, 
and the other subscripts k = R and H denote the observed pass- 
bands.' The ratio between the observed fluxes measured in 
the individual passbands are 

Fr _ AFr s + Frx Fr o _ Fr s+Frx 

Fh~ 



(2) 



AFh,5 + Fh;l Fh.q Fh.s + ,l 
where F/^s and fj. l represent the fluxes from the source star 
and lens, respectively. With some algebra, the color shift is 
represented as 

'(A + r]R)/(A + r]H)' 



A(R-H) = -2.5\og 



(3) 



(l + riR)/(l + r]H) 
where ijic denotes the ratio between the fluxes from the lens 
and source, i.e. 



m ■ 



(4) 



Fh,l 

Fr,s ' ■Ph,s 
The sign of the color shift is such that it becomes negative 
when the apparent source star becomes bluer during lensing 
magnification and vice versa. 

In the limiting case of a very low magnification (A 1.0), 
the amount of the color shift is approximated as 



A(/?-//)~-2.51og 



1- 



(5) 



(1 + 77r)(1+77h). 

where e = A - 1 .0. As A ^ 1 .0, e ^ 0, and thus A{R - 
H) 0, implying that the color shift expected from a low- 
magnification event is small. 

In the opposite limiting case of a very high-magnification 
event (A ^ 1.0), the numerator on the right side of equa- 
tion (O become (A + tir) / (A + rin) 1, and the color shift is 
approximated as 



A(/;-//)-2.51og 



(6) 



If the lens is very faint (i.e., rjk <C 1.0), equation Q is further 
simplified as 

A(R-H)^2.5log[l + (7^R-7^H)]- (7) 
Considering that the lens/source flux ratio is small, rjR-rjH 
1.0 and thus the expected color shift is also small. This im- 
plies that high magnification is not a sufficient condition to 
produce noticeable color shift. To produce noticeable color 
shift, it is additionally required that the lens should be bright. 
In the limiting case where the lens dominates over the source 
star in brightness (i.e., rjk'^^ -0), the color shift is represented 
as 

A{R-H)^2.5\og('^\={R-H)s-(R-H\, (8) 

where {R - H\ and {R - //)§ represent the colors of the lens 
and source, respectively. This implies that the color shift cor- 
responds to the difference between the colors of the lens and 
source. Under this condition, the color shift is maximized. 

' We choose R and H passbands because the MPF survey plans to observe 
in these bands. See details in § 4. 



3. SIMULATION 

We estimate the color-shift distribution of events ex- 
pected from future space-based lensing surveys by conduction 
Bayesian simulation of galactic microlensing events under the 
observational conditions of the MPF survey. For this simula- 
tion, it is required to model of the physical parameters includ- 
ing the mass density and velocity distributions of matter along 
the line of sight toward the field, lens mass function, and lu- 
minosity function of source stars and lenses. We conduct the 
simulation according to the following procedure. 

The survey is assumed to be conducted toward the galactic 
bulge field and the locations of the source and lens along the 
line of sight are allocated based on the spatial mass density 
distribution. Following the specification of the MPF survey, 
we set the center of the survey field at {l,b) = (l°.2,-2°.4). The 
spatial mass density distribution is modeled by adopting that 
of Han & Gould (2003). In the model, the density distribution 
of bulge matter is based on the deprojected infrared light den- 
sity profile of Dwek et al. ( 1995), specifically model G2 with 
^max = 5 kpc from their Table 2. The disk density distribution 
is modeled by a double-exponential disk, i.e. 



p(^,z) = Poexp 



R-Rq ^ \z\ 
Hr k 



(9) 



where {R,z) are the galactocentric cylindrical coordinates, 
R{) = 8 kpc is the distance of the Sun from the galactic cen- 
ter, pq = 0.06 Mq pc"-' is the local mass density, and Iir = 3.5 
kpc and h~ = 325 pc are the radial and vertical scale heights, 
respectively. 

Once the locations of the lens and source star are set, we 
then assign their velocities based on the velocity distribution 
model. The model velocity distribution has a gaussian form. 
The mean and dispersion of the bulge velocity distribution 
are deduced from the tensor virial theorem. The disk ve- 
locity distribution is assumed to have a flat rotation speed of 
Vc = 220 km s"' and velocity dispersions along and normal to 
the disk plane of a\\ = 30 km s~' and a±^ = 20 km s"', respec- 
tively. The resulting distribution of the lens-source transverse 
velocity from the combinations of the bulge and disk veloc- 
ities is listed in Table 1 of Han & Gould (1995), specifically 
the non-rotating barred bulge model. 

We assig n the source star brightn ess by using the luminosity 
function of lHoltzman et alJ (Il998l) . The Holtzman luminosity 
function is expressed in terms of V band. Once the V-band 
absolute magnitude of each star is chosen from the luminos- 
ity function, we determine the R an d // band m agnitude by 
using the color information listed in ICoxl (1 1999 1). The appar- 
ent magnitudes are then determined considering the distance 
to the source star and extinction. The extinction is determined 
such that the source star flux decreases exponentially with the 
increase of the dust column density. The dust column density 
is computed on the basis of an exponential dust distribution 
model with a scale height of = 100 pc, i.e. oc tx^{-\z\/h^) 
(iDrim mel & Sperg"elll2001t) . We set the normalization of the 
extinction in the individual passbands such that A^ = 1.8 and 
Ah = 0.8 for a star located at a distance Ds = 8 kpc toward the 
target field. 

The mass of t he lens is assigned based on the mass function 
of lGouldl(l2000l) . The model mass function is composed of not 
only stars but also brown dwarfs (BDs) and stellar remnants 
of white dwarfs (WDs), neutron stars (NSs), and black holes 
(BHs). The model is constructed under the assumption that 
bulge stars formed initially according to a double power-law 
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distribution of the form 
dN 
dM 



= k 



M 



0.7 Mq 



(10) 



where 7 = -2.0 for M > 0.7 Mq and 7 = -1.3 for M < .7 Mq. 
These slopes are consistent with lZoccalli et al.1 (120001) except 
that the profile is extended below their lower limit of 0. 15 Mq 
down to a BD cutoff of M = 0.03 Mq. Based on this initial 
mass function, remnants are modeled by assuming that the 
stars with initial masses 1 Mq < M < 8 Mq, 8 Mq < M < 
40 Mq, and M > 40 Mq have evolved into WDs (with a mean 
mass (M) ~ 0.6 Mq), NSs (with (M) - 1.35 Mq), and BHs 
(with (M) ~ 5 Mq), respectively. Then, the resulting mass 
fractions of the individual lens populations are stars : BD : 
WD : NS : BH = 62 : 7 : 22 : 6 : 3. Once the mass of the lens is 
set, its brightnesses in R and H bands are determined by using 
the mass-M y relation and the color information presented in 
ICoxl ( Il999l) . We note that BD and remnant lenses are dark and 
thus events caused by them do not produce color shift. 

With all the assigned physical parameters of the lens and 
source, we then produce lensing events. We do this under 
the assumption that the impact parameters of the lens-source 
encounter are randomly distributed. Since the lens flux does 
not participate in lensing magnification, it lowers the apparent 
magnification of the event, i.e. 



A = 



Mo + 2 



M0(M0 + 4)1/2' 



(11) 



where mo is the impact parameter of the lens-source encounter 
normalized by the Einstein radius of the lens. We define de- 
tectable events as those with apparent magnifications higher 
than a threshold value of A± = 3/5 that corresponds to the 
magnification when the source star just enters the Einstein 
ring of the lens, i.e. mo = 1.0. The apparent magnification 
varies depending on the observed passband. We use the ap- 
parent magnification measured in J band as the criterion for 
event detectability because it is the main passband of the MPF 
observation (see below). We assume that there is no blended 
flux other than the flux from the lens. We restrict that the 
survey monitors stars with apparent magnitudes brighter than 
J = 22.5, which corresponds to early M-type main-sequence 
stars. 

Once events are produced, the rate of the individual events 
are weighted by the factor 



r cx p(Ds)Dip(DOrEV, 



(12) 



where p(D) is the matter density, Dl and Ds are the distances 
to the lens and source, respectively, the factor is included 
to account for the increase of the number of source stars with 
the increase of distance, is the Einstein radius which corre- 
sponds to the cross-section of the lens-source encounter, and 
V is the lens-source transverse speed. The Einstein radius is 
related to the physical parameters of the lens and source by 



rE = 



4GMDdDs-DL) 



Ds 



(13) 



From the simulation, we find that 66% and 34% of the total 
events are caused by the lenses located in the bulge and disk, 
respectively. We also find that events caused by stellar lenses 
comprise 48% and others are produced by dark lens compo- 
nents. 

The MPF survey plans to observe in the wavelength range 
~ 600 - 1700 nm in three passbands. These consist of a very 



wide passband spanning the full 600 - 1700 nm range and the 
other two narrow bands at each end of the range (iBennett et alJ 
I20P). The very wide filter has no equivalence in standard fil- 
ter, but we assume it corresponds to J filter because the cen- 
tral wavelength matches that of the J filter The two narrow 
filters are approximated by the R and H filters. The observ- 
ing strategy for MPF is to do most of the observations with 
the very wide filter and a small number of observations in the 
R and H filters. We choose color is measured on the images 
taken in R and H filters because larger color shift is expected 
as the wavelength gap between the filters becomes wider. We 
assume that MPF can detect 14.3 and 5.0 photons s ~' for a 
22 ma gnitude star in the R and H bands, respectively (iBessell 
1979; Campins et al.]fT985h . 

For the individual events, we calculate the maximum color 
shift by using equation (O. The MPF plans to take images 
with 100 sec exposure at a 2 min interval and combine the 
images to make a 10 min exposure image, resulting in a data 
acquisition rate of ~ 3-5 images hr~' or ~ 100 images day"' . 
If 20% of the observations are done in the narrow filters, each 
of the R and //-band image is obtained every 2.4 hrs. This 
rate is frequent enough to cover the peak of nearly all lensing 
events during which the color shift is maximized. We assume 
that color shift is detectable only if the color shift measured at 
the peak of the event is five times greater than the uncertainty 
in the color-shift measurement, i.e. 5a detection. 

We determine the uncertainty of the color shift measure- 
ment as follows. Let us define TZf as the ratio between the 
flux ratios of an event measured in two passbands during and 
before lensing, i.e.. 



Tip 



Fr/Fh 
FRfi/Ffffl 



(14) 



According to the standard error analysis, the fractional uncer- 
tainty in the measurement of TZ is estimated as 



Tip 



Fr/Fh 



^Frq/Fh^O 

FRfi/Fnfi 



1/2 



where the individual terms on the right side are 
\Fr/Fh J \Fr J \Fh 



(15) 



(16) 



and 



'^Fr^ifhs V ^ V "^^" o V I f \ ^ (17) 
Fr.{)/Fh.()J \Frx) J \FHfl J 



Under the assumption that the photometry is photon domi- 
nated, i.e. ap =F^I'^, equations (fT6b and ( fTTb are approxi- 
mated as 



f ^ Fr/Fh 

\ Fr/Fh J Fr Fh FrFh 



1^1 Fr + Fh 



(18) 



and 



^FRa/F,,o 



1 I ^1 _ Frj+Fh.o ^ ^jg^ 



^FRfi/Ff].() J Fro F[f,Q pRoFufi 
Then, from equation (fTSb . ( fTSl ). and iT% , one finds 



Fr + Fh ^ Fra + Fh,o 



TZf \ FrFh FraFha 



1/2 



(20) 
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Fig . 1 . — Color shift distribution of events expected to be detected from the 
future MPF lensing survey. Darker shade region represents the distribution 
for events where the color shift can be detected at 5cr level. 

The uncertainty in units of magnitude is related to the frac- 
tional uncertainty cttj^ /TZf by 



crA(R-H) = 2.5loge 



TZf 



1.1 



TZf 



Then, the color-shift uncertainty is expressed as 



CrA{R-H) ~ 1.1 



Fr + Fh ^ Fr o + Fh,q 
FrFh FnfipHfi 



4. RESULT 



1/2 



(21) 



(22) 



In Figure[Tl we present the color shift distribution of events 
expected to be detected from the MPF survey. From the distri- 
bution, we find that color shift can be measured for 12% of 
the total events. Considering that events produced by stellar 
lenses comprise ~ 48% of the total event, this corresponds to 

23% of all events produced by stellar lenses. We find that 
the fractions of color-shifted events among the disk and bulge 
events are ^ 18% and 9%, respectively. The fraction for the 
disk-lens events is higher because the lenses are located closer 
to the observer and thus they are brighter 

We investigate the properties of the color-shifted events. 
First, we examine the brightness of the source stars and lenses. 
In Figure |2] we present the //-band brightness distributions 
of the sourceH-lens combined images (top panel), and separate 
source stars and lenses (second and third panels, respectively). 
From the distributions, we find that the color-shifted events 
have apparent brightness distributed in the range of 17 < // < 
21. We also find that the source and lens stars of color-shifted 
events have similar brightness ranges of 18 < // < 22. This 
implies that color shift occurs when the lens brightness is 
equivalent to that of the source star. This can also be seen 
in the bottom panel of Figure |2] where we present the distri- 
bution of the //-band lens light fraction. We find that ^45% 
of the color-shifted events will be produced by lenses brighter 
than source stars in the H band and ~ 1 1 % of them will have 
lens light fractions greater than Fh.l/(Fh,s+Fh,l) = 0.9. 

Second, we investigate the magnifications of the color- 
shifted events. For this, we construct the distribution of the 
impact parameters of the lens-source encounter. The distribu- 
tion is presented in the top panel of Figure [3] From the distri- 
bution, we find, as expected, that color-shifted events tend to 
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Fig. 2. — Top panel shows the distribution of the apparent brightness of the 
combined image of the source star and lens. Second and third panels are the 
brightness distributions of the separate source stars and lenses, respectively. 
Bottom panel shows the distribution of the //-band lens light fraction. The 
lens light fraction is 1.0 when the observed flux is purely from the lens. Dark 
shade regions represent the distributions for events where the color shift is 
detectable at 5(t level. 

have high magnifications. We estimate that the fraction of the 
color-shifted events with impact parameters less than mq = 0.1 
(or A > 10) is ^31%. On the other hand, the fraction of the 
color-shifted events among those with mq > 0.5 (or A < 2.4) 
is merely ^ 10%. These results imply that color shifts will be 
detected mainly from high-magnification events. 

Third, we investigate the event time scales. The event time 
scale is defined as the time required for the source star to cross 
the Einstein radius of the lens, i.e. fe = ^e/ v. The middle panel 
of Figure |3] shows the distribution of the time scale. By com- 
paring the distribution with that of all events, we find that the 
time scales of the color-shifted events tend to be longer than 
the events without color shifts. The most frequent value of 
the time scale for the color-shifted events is ^ 18 days. On 
the other hand, the mode value of all events is '--^ 6 days. We 
find this tendency is because the lenses of the color-shifted 
events tend to be brighter and thus heavier than the lenses of 
events without color shifts. We find that the lenses responsible 
for the color-shifted events have masses located in the range 
of 0.2 Mq < M < 1.0 Mq. These correspond to stars with 
spectral types from G down to mid M-type main-sequence 
stars. Color shift measurement for events produced by late 
spectral-type lens stars is possible because the MPF survey 
will observe in near infrared through which the lens/source 
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Fig. 3. — Distributions of tlie normalized impact parameters of the lens- 
source encounter (top), event time-scale (middle), and lens masses (bottom) 
of events to be detected by the future MPF lensing survey. Dark shade region 
represents the distribution for events where the color shift is detectable at 5cr 
level. 

flux ratio is not negligible. 

The selection effect of color-shift measurements that favors 
events with higher magnifications and longer time scales im- 
plies that the fraction of planetary microlensing events with 
color-shift measurements would be substantially higher than 
the average value of all events. This is because planet detec- 
tion has a similar selection effect of preferring events with 
higher magnifications and longer time scales (.Han .2007 ). 



Therefore, it would be possible to characterize the host stars 
of planets from color-shift measurements for an important 
fraction of planetary events. For some cases, the source and 
lens stars may have similar colors resulting in small color 
shifts. In this case, the lens stars can still be detectable with 
space-based fo llow-up observations condu cted several years 
after the event dBennettllAnderson & Gaudil) . 

5. CONCLUSION 

We investigated the usefulness of future space-based lens- 
ing surveys in measuring color shifts of microlening events 
caused by the flux from the lens. By conducting simulation 
of galactic microlensing events based on the specification of 
the proposed MPF survey, we estimated that the shift in the 
color of R-H would be measured at 5a level for ^ 12% of 
events that would occur on source stars with apparent mag- 
nitudes brighter than J = 22.5. Color-shifted events tend to 
have high magnifications and the lenses would have bright- 
nesses equivalent to those of source stars. We estimated that 
^ 3 1 % of the color-shifted events would have magnifications 
A > 10 and the lenses of ~ 45% of the color-shifted events 
would be brighter than source stars. The time scales of the 
color-shifted events would be longer than those without color 
shifts because these events tend to be brighter and thus heav- 
ier than those without color shifts. We estimated that the most 
frequent time scale of the color-shifted events would be ^ 18 
days, while the mode value of all events would be ^ 6 days. 
The spectral type of the lens stars of the color-shifted events 
would extend down to mid M type because the MPF survey 
will observe in near infrared through which the lens/source 
flux ratio is not negligible. Being able to extract information 
about the lens from the measured color shift, future space- 
based lensing surveys would be able to better constrain the 
lens population of galactic microlensing events. 

This work was supported by the Basic Research Fund of 
Chungbuk National University. 
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